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Abstract 

We extend our previous determination of the thermodynamic pressure of the Standard Model 
so that the result can be applied down to temperatures corresponding to the electroweak crossover. 
This requires a further resummation which can be cleanly organised within the effective theory 
framework. The result allows for a precise determination of the expansion rate of the Universe 
for temperatures around the electroweak crossover. 
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1 Introduction 


In a recent paper [1] we calculated the pressure of the standard model at high temperatures to three 
loops, or to order in the coupling constants. That work followed a long series of calculations 
dedicated to understanding the perturbative expansion of the pressure of gauge field theories at high 
temperatures. Especially, within QCD such calculations have been important: the computation of 
the coefficients of the expansion in g at high temperature has a long history and the result is known 
today up to the last perturbatively calculable term of order g^ In g [2-4], marking an endpoint to 
an impressive set of computations carried out in [5-10]. An optimal approach to calculating those 
coefficients is the effective theory method [11], based on asymptotic freedom and on separating the 
relevant mass scales: ttT, the electric scale toe = gT and the magnetic scale uim = g^T. 

An obvious drawback of those calculations within QCD is that they cannot be extended to study the 
QCD phase transition since the coupling grows large and thus any perturbative calculations become 
unreliable. That is not the case for the electroweak sector of the standard model. The Landau 
pole related to the weak interactions corresponds to a length scale 1/Aew ~ 10® m and therefore 
the confining effects can be expected to be negligible. As a result, calculating the properties of 
the electroweak phase transition using, for example, perturbative 1-loop [12-14] and 2-loop [15-17] 
effective potential calculations was possible for small Higgs masses. A better approach, well defined 
for large Higgs masses as well, is to perturbatively match the full 4-dimensional theory to an effective 
3-dimensional theory [18] and then to numerically solve the phase diagram from the effective theory 
using lattice Monte Carlo techniques [19-23]. Such studies show that the phase diagram has a first 
order line which ends in a 2nd order critical point of Ising universality class [24]. For experimentally 
allowed Higgs masses the transition is a crossover. Similar techniques have been used to solve the 
phase diagram also when the external U(l) magnetic field [25] or the chemical potentials related to 
the baryon and lepton numbers [26] are nonzero. Solving the phase diagram with numerical studies 
of the full 4-dimensional theory has also been achieved [27]. 

Although the properties of the electroweak phase transition are known today, computation of the 
pressure of the theory has been missing. However, the ongoing and future measurements of cosmic 
microwave background radiation (WMAP and Planck) allow for a precise determination of, for exam¬ 
ple, the WIMP relic density, which depends on the equation of state of matter in the early universe. 
A precise determination of the pressure is needed to match the accuracy of the observations [28]. 
In [1] we performed this calculation when the temperature of the system is high. In the present paper 
we will extend that result to temperatures close to the critical temperature of the electroweak phase 
transition. Critical temperature and the phase transition should be understood in perturbation theory 
framework, where there is always a first order phase transition. We will use these terms throughout 
this paper, although the actual transition is just a crossover. 

The present computation requires a reorganization of the effective theories. Close to the transition 
the fundamental scalar (which drives the transition) becomes light with respect to the electric scale 
and thus in a consistent calculation we have to formulate an additional effective theory by integrating 
out the adjoint scalars. The remaining theory contains just the fundamental scalar and the gauge 
fields. The pressure will then be composed of three parts: the contribution from fermions and the 
nonzero Matsubara modes of bosons ~ 1 + g"^ + g*{lj€ + 1), from the adjoint scalars Aq 

and Bq, pmi/T ~ + g\m^{\lt -|- 1) -I- g^mYy{l/e -I- 1) and finally from the fundamental scalar 

PM 2 /T ~ TOg -I- g\m^(\lt + 1) where ~ g^T’^ and m§ ~ At each step of the calculation 

there will be 1/e poles from ultraviolet and infrared divergences which will not cancel until all the 
contributions are summed together. 

The paper is organized as follows: in section |21 we explicitly define the theory we are working with, 
fix various conventions and briefly review the method of dimensional reduction as applied to the 
present case. Sections 13 and 01 contain the essential calculations and in section[3we discuss the result. 
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2 Basic setting 

We consider the SU(3)c x SU(2)l x U(l)y standard model with np = 3 families of fermions and 
ns = 1 fundamental scalar doublets, and evaluate the pressure of this theory at temperatures slightly 
above the electroweak phase transition. The theory is specified by the Euclidean action (in the units 
h = c=l) 


r/3 


5 = 


dr / d^xC 


( 1 ) 


£ = 


where 


+IlPIl + &Rpe.R + qLpqL + urPur + dnpdR + igy {qLT^^*tR — , ( 2 ) 


G% = d^A:-d.Al + ge‘^’^^AlAl, 


F„, 


= df,By - d„B^, 




Ph 

P^R 

pqh 

pUR 

pda 


5^$ _ 




= 7m 


dPL - + \bPl 


= 7 m+ ig'B^eR), 


13 


*5 


= 7m ^M'l'i - - -^B^qL - igsGp'°-qL 


= 7m 


= 7m 


6 

2 z (7^ 

d^UR - ^Bfj,UR - igsCpr°‘UR 

df^dR + '-^B^dR - ig.Cpr^^dR 


( 3 ) 


Here 5^ and (7“ are gauge bosons of weak-, hyper- and strong interactions, respectively; $ is 
the fundamental scalar doublet; p and cr are the left-handed lepton doublets and the right-handed 
lepton singlets (wrt. weak charge), and gz,, ur and dR are the left-handed quark doublets and the 
right-handed up and down -type quark singlets. The Yukawa coupling is taken into account for the 
top quark only. Summation over different families is assumed. Also, d = 3 — 2e in dimensional 
regularization, v = 0 ,..., d. The gamma matrices are defined in Euclidean space so that { 7 ^, 7 ,^} = 


{75,7m} = 0 and Tr757^7,,7p7,^ oc 


The color indices are a = l,...,dA for the weak 


interaction and a = 1, ■■■, — 1 for the strong interaction. The different group theory factors for 

SU(A^) with generators T°‘ are defined as: 
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For SU(2) with T“ = r“/2 they are Tp = 1/2, Cp = 3/4, Ca = 2, dA = 3 and dp = 2. 

The theory in Eq. @ contains six couplings that run with the renormalization scale: gauge couplings 
g', g and the fundamental scalar quartic self-coupling A and its mass parameter and gy- We fix 
the values of these couplings at the scale /r = mz according to their tree-level relation with different 
physical parameters as in [1], where the 1-loop running of the parameters is also given. We employ 


a power counting rule A 
1/2 ^ g^T\ 


9‘ 


98 
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g'^ and assume the temperature to be such that 
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Figure 1: The ratio of the fundamental scalar thermal mass to the SU(2) adjoint scalar (Debye) mass. 
The expressions for the masses are given in Appendix m and the regularization scale is chosen as 


A = 2ttT. 


The physical observable we are studying is the pressure, defined by 

p{T) = lim —In f'DA'Dip'Dtp'D^ exp {—S). ( 6 ) 

V —>00 V J 

It is normalized such that the (real part of the) pressure of the symmetric phase vanishes at T = 0. 
Other interesting variables, such as entropy and energy densities s(T) and e{T), can then be evaluated 
using standard thermodynamic relations, s{T) = dp/dT, e(T) = Ts{T) — p{T). The purpose is to 
calculate the pressure up to, and including, order + lng)T^, employing the power counting rules 
above. For high temperatures, this calculation was performed in [1]. However, close to 

the phase transition, which is often the most interesting region, our earlier computation is not valid, 
since the (thermal) fundamental scalar mass is much smaller than the Debye masses and introduces 
another hierarchy of scales which needs to be sorted out. The ratio of the two mass scales at different 
temperatures is shown in Fig. ^ which illustrates that there is a range of temperatures in which 
m 3 mD. The expression used for m\ in Fig. ^ is Eq. (I4till . renormalized in the MS scheme and e 
set to 0 . 

It is well known that simply evaluating all 1-, 2- and 3-loop vacuum diagrams in perturbation theory 
fails because of infrared divergences. Instead, one has to separate the contributions of different scales 
into successive effective theories [11], where all the large scales are integrated out one by one. First, 
we integrate 

p{T) = pEiT) + ^ In JvAkVAoV^ exp (-,Sei) , (7) 

where S'ei contains only the static Matsubara modes of the gauge bosons and of the fundamental 
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scalar (Higgs) field. The contributions of the nonzero Matsubara modes and fermions to the pressure 
show up as the matching constant pe (App. © and in the parameters of Sei (App. 10 ). The spatial 
(magnetic) gauge field components remain massless, while the temporal component gets a thermal 
mass Tod ~ gT. The theory defined by S'ei can then be viewed as a 3d gauge theory with adjoint and 
fundamental scalar fields. 

The effective theory thus obtained still contains contributions from several scales, TOq ~ g'^T^, 
m\ and the magnetic scale g'^T. We restrict ourselves to temperatures close to the phase transition, 
where m| is small compared to the Debye masses and assume m\ < g^T^, motivated by Eq. 11411 . We 
integrate out the scale gT, i.e. the fields Aq and Bq, and are left with 

p(r) =pe(T)+pmi(T) + ^In JvAkV^exp{-SE2) ■ (8) 

This is a 3d gauge theory with a fundamental scalar field. The largest mass scale in this theory is 
the fundamental scalar mass < g^T^, so this theory contributes at ~ g^^^ or higher order, and 
we need to evaluate only 1- and 2-loop vacuum diagrams from this theory. The potential infrared 
problems can be isolated to another effective theory by integrating out the fundamental scalar, 

p{T)=pE{T)+pMiiT)+pM2{T)+'^\n J'DAk'D^exp{-SM) ■ (9) 

The remaining effective theory contains only the (massless) spatial gauge fields. Therefore the only 
mass scale of the theory is provided by the 3d gauge coupling and is of the order g^T and consequently 
the contribution of this theory to the pressure is of the order g^. The final result of our calculation 
can then be written as 

p[T) = pe(T) +pmi(T) +pm 2 (T) +pQCD(r) + 0{g^-^T^), (10) 

where pqcd can be taken from [2,8-10]. One-loop quark diagrams are included in pE, so they must 
be subtracted from pqcd- The pressure pe was computed in [ 1 ], and is given in Appendix 1X1 whereas 
Pmi and pm 2 are computed in Sections 0 and 0 respectively. 


3 Calculation of the pressure pMi 


In this section we integrate over the scale gT. This means that we need to integrate out the adjoint 
scalars Aq and Bq, leaving an effective theory with only gauge bosons and the fundamental scalar 
field. 

The 3d theory of gauge fields with fundamental and adjoint scalars is given by 


5ei = J d^* + (A$)^(A$) + + A3($t$)2 

+ 2(AA[))^ + + '^^a{AqAq)'^ + -{diBo}'^ + -iji'^BqBq 

+h3$^$AgAg + K'q^^^BqBq - i53ff3^o$^Agr“$ 


( 11 ) 


where the parameters were computed in [1] to required order and are listed in Appendix 151 

To find out the contribution of the adjoint scalars to the pressure we need to calculate all the 
vacuum diagrams in the theory m containing these fields. Since both Aq and Bq are massive, there 
are no infrared divergences related to this integration. The mass of the fundamental scalar m 3 is 
parametrically smaller (m| ^ g^T^) so we treat it as a perturbation and expand the integrands in 
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Figure 2: Diagrams contributing to pmi- The dashed lines correspond to the fundamental scalar, the 
dot-dashed lines to the adjoint scalars, the wavy lines to the gauge fields and the dotted lines to the 
ghosts. 


toI/toq. In fact, the leading m 3 = 0 order is sufficient to our computations. The required diagrams 
are shown in Fig. Inland are evaluated in [1]. 

The result, with np = 3 and ns = 1, reads 
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The divergences and the scale dependence in 0{g‘^) cancel against those in pu, and the ones in 0{g^) 
against the sunset diagram in pm 2 , Eq. 113. 


4 Calculation of the pressure pm 2 

The theory containing just the gauge fields and the fundamental scalar is defined by 

5e2 = y'd3a;|iG“G“ +iF,,F,,-P(A‘i>)^(A$)+m2<i>t$ + A3($t$)2|^ (13) 

where the field strength tensors G“ and Fij and the covariant derivative are defined as before 
(with couplings 53 and ^ 3 ). The integration from Sei to S'e 2 will not introduce any new corrections to 
the couplings to get the pressure to the required order and therefore (gf, g'^, A 3 ) = (g|, g'^, A 3 ). The 
scalar mass does, however, get an additional contribution (calculated in [18] apart from the term 0 (e) 
which is needed in the present calculation) coming from the tadpole diagrams with adjoint scalars at 
the loops: 

dAhzmu + ^g'im'u 
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Figure 3: The diagrams required to calculate pm 2 - The dashed lines correspond to the fundamental 
scalar and the wavy lines to the gauge fields. 


1 

27r 


dAhsUlB 


+ In 


\ 

2mD J 


12 / 
753 


1 + ln 


2 m'r, 


Oig% 


(14) 


Since the pressure is to leading order given by pm 2 /T ^ m 3 and since we assume the temperature 
to be such that rhg ~ g^T^, we see that corrections of order g'^ to the scalar mass would contribute 
to pressure at order g^'^ and we will neglect those. Similarly, we can neglect the g'^ corrections to 
m§. The running of the mass m§ (and hence of m§ as well), which sets in at the order g'^ does not 
therefore influence the present calculation as it did the high temperature calculation. Furthermore, 
even though m§ is assumed to be of the order g^ near the phase transition, the expression for m§ 
does contain individual terms of the order g^ which, when summed up, effectively cancel to order 
g^ (assuming proper temperature region). The individual terms of order g^ must, however, be taken 
carefully into account when cancelling the 1 /e poles in the pressure. 

The pressure is obtained by calculating the diagrams in Fig. |21 Three-loop diagrams would con¬ 
tribute to order ® and are neglected now. The result is 


PM2 

~T~ 



m| 

(47r)2 


dp{dF + 1)A3 -I- -dp 



-h3-f 41n 



The 1/e poles cancel against those coming from pp and pmi- 


(15) 


5 Numerical results 

In this section we will consider some numerical implications of the result, given by Eq. HI OH . The 
convergence and the scale dependence of the perturbative expansion were already studied in [ 1 ] so 
now we will just concentrate on the difference between the results of that paper and the present 
paper, which are significant when the fundamental scalar becomes light. The Higgs mass is taken to 
be 130 GeV and the W mass 80 GeV. 

It is instructive to first look at a simpler SU(2) -|- fundamental Higgs model for which the per¬ 
turbative expansion is better behaved. In Fig. 01 we have plotted the pressure of this theory in the 
temperature region where the fundamental scalar becomes much lighter than the adjoint scalars (see 
Fig.nj using both the calculation valid at high temperatures (pht) and the calculation valid near the 
phase transition (ppx), normalized to po = (ideal gas pressure of SU(2) gauge bosons and a 

massless fundamental scalar). It can be immediately noticed that there is a temperature region where 
the high temperature calculation is not well behaved but becomes singular. This is due to the terms 
of the type g'^ra'^/m^ in the expansion of pht which become singular when m 3 —>0. In the present 
calculation which takes into account that <C such terms will not appear and therefore ppT is 
seen to be well behaved.^ It is worth noting, however, that the temperature region where the singular 
behavior is manifest is very narrow. That is due to the small gauge coupling: the singular terms are of 

®Note that as frig becomes negative, ppx will develops an imaginary part since the symmetric phase, where the 
calculation is performed, is not stable anymore. The imaginary part can then be related to the rate of decay of the 
symmetric phase [29]. In the figures we have plotted the real part of the pressure. 
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Figure 4: The pressure of the SU(2) + Higgs model, ppt is the computation which takes into account 
that TO§ <C pht is the high temperature result. 

the order and their contribution to the pressure vanishes fast as one moves away from the singular 
point. More specifically, the singular terms behave as: 

^singular _ 135 g^ ^ j Tq H fii 

PO " 4096^5 ^1^2 

where Tq is the location of the singularity and ST = T — Tq. Since the curve behaves as ^ l/V^ 
near the singularity and since the prefactor is so small, the contribution from the singularity becomes 
small very fast as ST is increased. 

If one excludes this narrow region around the singularity, the pressure ppT is smaller than the pres¬ 
sure pht roughly by constant xT^. The relative difference between pht and ppx is small, pht/ppt — 
1 ~ 0.02 but it nevertheless is of the same order of magnitude as differences in the pressure at dif¬ 
ferent orders in perturbation theory and hence it is important to take this effect into account when 
calculating the pressure to high accuracy. This is related to the reorganization of the effective theo¬ 
ries. Keeping the fundamental scalar mass, as defined in Eq. m, unexpanded when inserted to the 
expression for the pressure pm 2 > Eq. m, effectively resums all the ring diagrams of the type shown 
in Fig. [3 to the pressure. The high temperature calculation, however, corresponds to expanding this 
mass in terms of h^mu/m^ in the expression for the pressure and to keeping just the lowest order 
terms of the expansion. At high temperatures this expansion can be made, since ~ g 

but close to the phase transition ~ 1 and the expansion is not possible. Moreover, even 

though differences in the pressure when calculated from the high temperature theory as compared to 
the correct calculation near the phase transition are small, this is not necesserily the case for other 
interesting variables such as the critical temperature and to obtain these correctly, one needs to use 
the theory that takes into account that the fundamental scalar becomes light. 
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Figure 5: The type of diagrams that are resummed to the pressure pm 2 - The dashed line corresponds 
to the fundamental scalar and the dot-dashed to the adjoint scalars. 


Although we do not study the specific features of the phase transition, which are already well 
established [18-27], we nevertheless plot the pressure of the broken phase as well in Fig. 01 to order 
to indicate the temperature where the phase transition takes place. In the previous sections the 
pressure was calculated in the symmetric phase and thus the result cannot be extended to the broken 
symmetry phase as such. However, we can make use of the effective potential calculations [15-17] 
and write pBp{(p) = psp — V{(p), where the effective potential is normalized so that H(0) = 0. The 
pressure in the broken phase to 0 {g^) is then given by 


ppp{TM = W-V + 


192 ‘ 


24 
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[mnipf + imcBipf + 9mw{pf) 


-b 
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6 mw(p)^ + 3 ( mwip) 
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+ 0{g\ 
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where mH{<p)'^ = 3X(p‘^ — and mw{p)^ = g^'pP 'are the zero tempera¬ 

ture masses of the particles, and g) = y>(r) (the expectation value of the Higgs field) is such that 
dppp!dipP' = 0. The critical temperature below which the pressure of the broken phase is bigger than 
the pressure of the symmetric phase, indicating that the symmetry is spontanously broken, is, for 
tuh = 130 GeV and mw = 80 GeV, about Tc = 215 GeV. 

In Fig. El we have plotted the pressures pht and ppx for the full standard model, normalized 
to po = 106.757r^T^/90. A similar structure can be seen to appear here: the high temperature 
calculation of the pressure developes a singularity when 0, but the temperature region where 

this is of importance is very narrow, while the calculation which takes into account that near the 
phase transition is seen to behave well. Excluding the singularity, the pressures are again 

observed to differ roughly by a constant xT"'’, but the difference is not in this case even as large as in 
the SU(2) + Higgs theory, especially when compared to the corrections introduced by each new order 
of perturbation theory. This reflects the fact that the fundamental scalar carries only a small part of 
the total number of degrees of freedom and hence changes in the scalar sector of the standard model 
do not lead to significant changes in the pressure. 


6 Conclusions 

In this paper we have extended the calculation of the pressure of the standard model at high tem¬ 
peratures carried out in [1] to the case when the temperature of the system is close to the critical 
temperature of the electroweak phase transition Tc. The previous calculation was not consistent in 
that region since it assumed that the fundamental scalar has a mass comparable to the masses of the 
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Figure 6: The pressure of the SM, ppx and pht are the present and the high temperature computation, 
respectively. 


adjoint scalars while, in fact, near the phase transition the fundamental scalar becomes light. This 
inconsistency manifests itself as an unphysical singularity in the pressure if one tries to apply the 
result of the previous calculation to temperatures near Tc- The calculation performed in this paper 
is done in a consistant manner, taking into account the lightness of the fundamental scalar near the 
phase transition and no such unphysical singularities remain in the final result. 

It is possible to go still a bit further in the perturbative expansion of the pressure, the next term 
would be of the order g^'^. Its calculation would require a two-loop matching of the fundamental 
scalar mass (order g'^) and a three-loop calculation of the pressure pm 2 , both of which are in principle 
manageable. Furthermore, one could determine the final perturbatively calculable term g^hig by 
calculating the vacuum energy densities of the three dimensional effective theories. However, since 
the convergence of the perturbative expansion is rather fast (as noted in [1]), especially for the SU(2) 
-I- Higgs theory, the calculation of those terms would not lead to large numerical differences in the 
pressure. 
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A The pressure 

For completeness, we will give here the result for the pressure pe as calculated in [1]. The general 
form for pE can be written as: 


Pe(T) = 
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- Ar , 1 19, A 

16 V e 4 AeT 

- Ar , 1 131 , A 

—Nr ( —I-In 


48 

1 


36 AeT 


^_13i,2 + ,+ I^ + i^ 

120 4 2 C(-l) 4 C(-3) 

6563 41, 23C'(-1) 

1080 20 ' 18 C(-l) 


49 C'(-3) 
^ C(-3) 


A 


—Nr In-In 2 4- 7 

6 V AeT ’ 

Cpdp 


12 


(iV.^-l)„,(|-l„2) 


1 

12 


[Nr - 1) np 


175 5, „ dp /35 , „ 

-ln2-4—-ln2 

288 9 36 132 


15 

~144 

dp 

—ns 
12 


{n; - 1 ) (in 


A__ii_ii,„2+2i:^-i:^ 

AeT 75 25 C(-l) C(-3) 


Cpdp /I A , 5 , , ,,C'(-1) 

e+^^"4^ + 3+9' + 2 ^ 


ns 


(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

( 31 ) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 
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C'(-I) 

C(-1) 


(39) 


OIEBv 

OLEXv 

aEYv 


OiEuv 


dp 

-3-"S 

‘'’■"=(‘"4^-1 + ^) 


1 o, A 5 



(40) 

(41) 

(42) 


B Parameters of S^i 


For the couplings, tree level values are sufficient for our purposes: 


gi 

= g^T, 

g'i 

= g'^T, 


A 3 

= at. 

\a 

= 0 {g^), 

(43) 

ha 

= h^T, 

h’a 

= y^T. 



For the mass parameters we need terms of order (\-\-e)-\- . These were computed in [1] and we 

just quote the results here. The adjoint scalars do not have divergences to this order, 


2 _ t~i2 

mo = T 


{Pei + PE2e + 0{^)) + ^ {Pe^ + 0{e)) + 0 (/) 


(47r) 


2 I Pe\^ + PEsgi + PEYgy 4 " Ps'g'^ + Pei>' 


, 2 \ 1 


2^2 


(44) 


Wd 


— T'^ I g'^ {P'ei + Pe 2 ^ + ^(e^)) + + ^(e)) + ^{g'^) 


i4 


g' 


(47r)2 


P'e\^ + Pssgs + Psygy + Psg + 


y 2 


(45) 


The fundamental scalar, on the other hand, has 1/e divergences at g"^ order. The bare mass reads 

T2 


mg (A) + Sm^ = 


4 I /4 I 4® 2 /2 ® \ 2 3 \ /2 I o\2 

xo I- Q —Q H” —Q Q — —H” 3A 

(47r)2e V 64^ 64^ 32^ ^ 4 ^ 4 ^ 


g^ g'^ 

1 + T ":—~PuA + -r. — ^Pi/B 


A 


;P. 


(47r)2'^""‘ ' (47r)2^"° ' (47r)2^‘^^ ' (47r)2 

+ [g'^iPAl + PA2p) + g'^iPBl + Pb2P) + A(/3ai + /?A2c) + gyiPYl + Py2() 

A 2 


gy 


■P, 




^4 /4 Xfp' \n'‘^ 

Paa + tj-^Pbb + tj-^Pab + tt-^Pax + tt-^Pbx 


(47r 


(47r)2 


(4^)2 


«2„2 n'2 2 2 2 \ 2 4 

I g gy p. , g gy ^ , 9 y „ 

T /. ,2 Pay T /. o Pby T /. _',9 Psv T /. _\9 pay T /. 9 Pw 


( 47 r) 


(47r)2 


(47r)2 


(47r)2' 

(47r)2 


(47r)2 


(47r)^ 


/3 a A 


(47r)2 


(46) 


where (5m| is the counter term in the MS scheme. 

The constants in the above expression, with the correct group theory factors substituted into them, 
are as follows: 


Pei = 2 + gHp + -ns 

R' 

Pei = g^F + gns 


(47) 

(48) 
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Pe2 


P'e2 

Pe3 


P'e3 


1 + 1 In — 

3 C(-l) 3 47rr 

/I iC'(-i) 1 , 

+ (6 + 3CFTy^3^^4^T 


2 4, „ 4C'(-1) 4, A \ np 

3 3 3 C(-l) 3 AttTJ 2 

A 


ns 


5 40 ..O , lO C(-l) 
9 9+9 C(-l) 


A 


10 1 


, 1 IC'(-I) 1, A \ 

' 4 + 2CFTy+2^^4^j 


ns 


20 88 88 , A 


1 32, „ 28 

-In 2 H- 

2 9 9 


+ 


4 16, „ 8 8 , A \ 


rip 


28, A 

— in- 

9 4eT 


+ — TT + 


11 


13 13, A 

—7 H-In- 

6 ' 6 AeT 


ns 


1 8 , „ 5 5, A 

-In 2 - 7 -In- 

9 9 9 ' 9 AeT 


100 
5 


400, „ 200 


npns 

A 


200 , 


TZp + 


5 1 1 1 ^ 

72 ~ Ts^' “ 18 4^ 


ns 


40, „ 25 25, A 


npns 


M 


np 


(49) 

(50) 


(51) 


(52) 


Pe\ = 

ns 

II 

— 2 np 

Pey = 

1 

~A 


1 1 

Pe' = 

--„p + -„s 

II 

2 ns 


Pe\ — 

ns 

(53) 

P'es = 

22 

(54) 

Pey = 

11 

~12 

(55) 

P'e = 

1 3 

-R + R 

(56) 

P'ev = 

2 ns 

(57) 


Pi^A 

Pv\ 

Pai 

Pbi 

P\1 

Pyi 


9, A 

— In- 

2 47r 


- 127 - 12 In 


A 


3 

16 

1 

16 

1 

2 

1 

4 


PuB 

PuY 

Pa2 

Pb2 

P\2 

Py2 


-In — 
2 47r 


— 12 In 2 — 67 — 6 In 


1 

4 ’ 
1 

12 
1 + 


3 C'(-1) 

8 C(-l) 
I C'(-I) 
8 C(-l) 
C'(-I) 


^In 


+ 


1 


A 

A 

A^ 
A 


■In- 


C(-i) 


+ ln 


8 AeT 
A 


47 T 


1 1 IC'(-I) 1 A 

2 " 2^^^+2 CRT ^ 2 ^^ 4 )^ 


(58) 

(59) 

(60) 
(61) 
(62) 
(63) 


Paa 


17 75 C'(-1) 5 5, A 

--—^ — — In-+ 

24 32 C(-l) 32 ^ 2 AttT 

12 16 C(-l) 4 ' 16 AeTJ 


+ 1 ^ 2 - 7 -In 
ns 



np 

T 


(64) 
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f3BB 


Pab 

Pa\ 

Pbx 

P\X 

Pay 

Pby 

PsY 

PXY 

Pyy 


9 C'(-l) 


+ 


22 C(-l) 

1 1 C(-l) 

36 16 C(-l) 


3 3, 

—7 H— In 
32 ' 8 


A 

4^ VlOS 


/5 5,„5 5, A\ 

I -1-In 2- 7 -In- ) 

18 18 ' 18 47rT J 


Up 


A 


- 7 -in 

12 ' 48 47rT 


ns 


, 15C'(-1) , 
16 C(-l) 
15_9C'(-1) 
4 


9 3, A 

16^+2^" 4^ 


A 


_9 

2 C(-l) 2 ^47rT 

A 


= 6 


5 3 C(-1) 3^^_ 

4 2 C(-l) 2 47rT 

6AA>-6. 


C(-l) 

3 3 , „ 

l6-8‘”^ + 

48 72 


17^ 


9, A 
- in ■ 


32 

-2 + — In 2 + 47 + 4 In 
o 


4ttT 

17, A 

— In- 

24 47 T 

A 


47 T 


= —3 In 2 —67 —6 In 


A 

4^ 


3, A 

7 H— In-. 

' 4 AttT 


(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 

(71) 

(72) 

(73) 

(74) 
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